NASA TECHNICAL NOTE 


VO 

«-o 

vO 

CO 







<C 



o ^ 


LOAN COPY: Rl 

AFWL [WISH 
KFRTLAND AF£ 


HIGH-ENERGY ION BEAMS USED TO 
ACCELERATE HYDROGEN PROPELLANT 
ALONG MAGNETIC TUBE OR FLUX 


bj Gerald W. Englert 

Lewis "Research Center 
Cleveland, Ohio 





NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


OCTOBER 1966 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 

iiiiiiiiiiinii 

□13QES5 

NASA TN D-3656 


HIGH-ENERGY ION BEAMS USED TO ACCELERATE HYDROGEN 

PROPELLANT ALONG MAGNETIC TUBE OF FLUX 

By Gerald W. Englert 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.00 



HIGH-ENERGY ION BEAMS USED TO ACCELERATE HYDROGEN PROPELLANT 

ALONG MAGNETIC TUBE OF FLUX 


by Gerald W. Englert 
Lewis Research Center 

SUMMARY 

An analytical study was made of the transfer of energy from a high -intensity beam of 
light-weight ions to a hydrogen target. The hydrogen was confined to a region bounded by 
a magnetic tube of flux. Under certain conditions the gradient of the magnetic field was 
quite effective in accelerating sizable mass flows of hydrogen which resulted in large in- 
creases of momentum over that of the impinging beam. Interaction between the various 
species of particles was studied in detail. 

The analysis was applied to the efflux of ions from a hypothetical 100-megawatt ther- 
monuclear reactor. The stream of ions escaping from the weaker mirror of a magnetic 
mirror system was represented by a high-flux-density monoenergetic unidirectional 
beam. Hydrogen was admitted into this beam slightly downstream of the weak mirror 
exit station. Thus, the hydrogen served as the propellant and the reactor as the energy 
source of a thermonuclear rocket. 

A range of magnet sizes and strengths, energy and cross-sectional areas of the 
beam, and hydrogen flow rates were studied. Energy transfer was most effective when 
the beam was concentrated on a dense closely confined target. Ionization and accelera- 
tion of the hydrogen was usually quite rapid near the station where the beam first im- 
pinged on the target. Beyond a distance of about one magnet diameter downstream of 
this station, the velocity was increased and density decreased to values that permitted 
little further energy transfer. The fractional ionization then stayed essentially constant; 
however, acceleration continued for another two magnet diameters because of the con- 
tinuing high gradient of the magnetic field that interacted with the remaining thermal 
energy in the propellant electrons. 

At low hydrogen flow rates only a small amount of energy could be transferred from 
the beam to the target. At high flow rates, loss was mainly due to power consumed in 
ionizing the hydrogen. The momentum (thrust) of the hydrogen could be increased to 
seven times that of the beam with an overall energy transfer efficiency of 25 percent. 



INTRODUCTION 


Numerous theoretical and experimental stopping -power studies (refs. 1 to 4) were 
made wherein an ion beam impinges upon and transfers its energy to a gaseous target. 

In the present study the target is confined to a slender region bounded by an axisymmetric 
tube of magnetic field lines. The flux tube is concentric with the axis of the magnetic 
field coil. Such a flux tube can form a type of "magnetic nozzle" as energy gained by 
the target gas is in part converted into a directed translational mode along diverging 
field lines. The gas initially requires a confining wall. Where required, the wall is 
congruent with the field lines. However, the gas that is ionized later can be magnetically 
confined, A net increase of momentum, or thrust, over that of the impinging beam alone 
could result as the additional mass of target material (propellant) is accelerated. This 
report investigates the conditions under which such results can be theoretically realized. 

Application is made to a conceptual thermonuclear rocket, the basic constituent of 
which is a reaction volume confined by a magnetic mirror system having one mirror 
slightly weaker than the other (fig. 1), A stream of high-energy helium 3 and deuterium 
ions is emitted from the weaker mirror (refs. 5 to 7) providing a small amount of thrust 
at a very high specific impulse. For interplanetary propulsion the thrust of such a 
system must be greatly increased by the addition of much more propellant (ref. 8). Hy- 
drogen is used as the additional propellant and is accelerated by interaction with the 
emerging reactor ions and the magnetic field of the local mirror . The hydrogen has 
previously been used to cool a high -temperature shield about the reactor (fig, 1) and is 
in a dissociated or monatomic state. The present analysis approximates the efflux of 
reactor particles from the weaker mirror by a monoenergetic unidirectional beam. 

(The mks system of units is used unless otherwise specified. ) 


/-Injection Of atomic / 
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ANALYSIS 


The analysis is based primarily on the theory of classical particles. In general the 
ensemble of particles considered consists of electrons, ions, and neutral atoms. Parti- 
cles initially in the beam include electrons plus one or more species of ions. The pro- 
pellant, or target, is initially in the form of monatomic hydrogen, some of which is soon 
ionized to form two additional species - electrons and protons. 

The behavior of the propellant species was studied as usual by conservation equa- 
tions of mass, momentum, and energy. Beam particles, however, were considered in- 
fluential in this study only as a source of energy. Preliminary calculations indicated that 
in the range of interest for propulsion the density and flow rates of propellant particles 
are much greater than those of the beam particles and that only the propellant species are 
of any consequence in the mass and momentum balances. The energy of the beam species 
and the density, velocity, and temperature of the propellant are determined as functions 
of distance along the magnetic field lines. The study is limited to a one -dimensional 
analysis in which the flow properties are uniform across any surface normal to the flow 
direction. 


Transfer of Energy from Beam Particles 

Beam composition . - The energy is assumed to be initially in the form of a high- 
intensity monoenergetic unidirectional beam of charged particles flowing along magnetic 
field lines. Such a beam undergoes a certain fractional neutralization because of charge 
exchange when it impinges on a neutral target gas (ref. 9). Very little momentum trans- 
fer occurs during a charge exphange event, and, thus, very little change in directions is 
expected for the incident particles. Thus, if the charged particles initially impinge upon 
the propellant as a unidirectional beam, they should remain so even if neutralized. Each 
atomic species of the beam is assumed to be monoenergetic at any given axial station in 
the magnetic flux tube. This assumption is consistent with the data of reference 10, 
which show that there is very little spread of energy about an average as a beam of ions 
transfers its energy to a plasma. 

Impingement of charged particles on a neutral target has been studied experimentally 
by many investigators. In these experiments the bombardment chamber was large or 
attached to an abundant reservoir of target material. The beam thus underwent a change 
in charge composition whereas the target material, for all practical purposes, remained 
the same. Herein, however, the target boundary is the same as that of the beam and also 
the flow of target material is limited. The target can thus undergo sizable changes in 
charge composition. 
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If the target material has a low density (within a certain rai^e) enough collisions 
occur to permit the beam to come to an equilibrium composition; yet the collisions are 
few enough that the energy of the beam is not appreciably degraded (refs. 4 and 9). After 
equilibrium is reached the beam composition is essentially independent of the distance 
that the beam travels. Experiments conducted under these conditions are useful for 
studying charge exchange and ionization cross sections as functions of beam energy. 
Herein, however, the target is often dense enough that the beam energy is constantly de- 
grading along its path of travel. Thus, both beam and target have equilibrium composi- 
tions, which, among other things, are functions of distance in the direction of beam 
travel. Consider, for example, a beam that is initially composed of doubly charged 
helium 3 ions, singly charged deuterium ions, and electrons (see APPLICATION section 
and appendix B). Shortly after impingement on an atomic hydrogen target the beam is 
composed of doubly and singly charged helium 3 ions, singly charged deuterium ions, 
helium 3 and deuterium neutrals, and electrons; whereas, the target is composed of pro- 
tons, electrons, and neutral hydrogen atoms. The fractions of the various species de- 
pend on the energies transferred and thus on the distance inside the target. A detailed 
study of this phenomenon would necessitate the inclusion of approximately 25 interactions 
in an already complex problem. 

The problem is made somewhat more manageable by restricting the analysis to 
beams which are composed predominantly of neutrals and singly charged ions. Hydrogen 
and its isotopes obviously fall into this category. Helium can also be included if the beam 
energy is kept below 100 keV because the probability of a doubly charged helium particle 
obtaining an electron from a target hydrogen atom far outweighs (by more than a factor 
of 10) any of the other charge -exchange probabilities for this energy range (ref. 9). 

Stopping power . - The energy exchange per unit distance per incident particle 
dEjj/dx is often written as the number density of the target material n^ times the total 
stopping power (l/n^)(dEj^/dx). The total stopping power can be expressed as the 
weighted sum of the partial stopping powers of each of the charge species plus the energy 
loss due to collisions involving charge exchange (refs. 3 and 4). Thus, for this restricted 
analysis 


. dEp 1 <^EpO . dEp+ 

-= = F — + F — + Loss due to charge exchange 

n^ dx R° n^ dx n^ dx 

where Fj^o and Fj^+ are the equilibrium fractions of neutral and singly charged ions 
in the beam, respectively. 

A conclusion reached in reference 1 is that the partial stopping powers for both 
neutral and singly charged hydrogen atoms are close to the same magnitude. With this 
premise and n^lecting energy loss due to charge exchai^e, total stopping power was 
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closely approximated by 


+F \ 1 _ 1 

dx \ H'*’ H°j dx dx 

This procedure also approximates reasonably well the total stopping power for helium in 
the energy range below 100 keV (see appendix C). 

Analytical stopping power expressions are reviewed in references 11 and 12. All 
the expressions are of the form 


"t 


dx 


» Constant) 



slowly varying function of 



where the subscripts t and R denote target particles and incident particles, respec- 
tively. The mj^/m^ term shows that the relatively heavy nuclei of the target particles 
receive very little energy in comparison with that received by the target electrons. 

The stopping power expressions also show that a beam of incident electrons having 
energies of the same magnitude as the incident ions would give a relatively insignificant 
amount of energy to the target. On the other hand, if the beam electrons were at an 
energy so much lower than that of the ions that the energy transfer rate from them was 
significant, then their energy content would be so low that only a negligible total amount 
of energy would be available from them. Thus, the electrons in the beam were assumed 
to serve only as producers of a neutral net flux of charge of the beam species (this will be 
discussed further in the APPLICATION section). 

Stopping power for species closely related to the one of interest herein has been de- 
termined experimentally (refs. 3 and 4). The theory of reference 2 gives satisfactory 
agreement with these measurements; therefore, it was applied to the singly charged ions 
impinging on neutral hydrogen atoms in appendix C. From equation (C3) 


^2 In ® when U < KEj^(l - W|,) 


dE 


R 


"n® 


4 


Wr 3 


'R' 


dx 


IBire^E 


R 


(‘ + *r) 


3/2 




In + 1 . 804 + ^ ^ 


w 


R 


V. 


Wjj 6 96 


when U > KEjj(l - Wj^) 


( 1 ) 


5 


where 




Y“e '^R 
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(me + mj^)2 


the quantities Ej^ and are the energy and mass of a single -incident or impinging- 
reactor particle, and n^^, Eg, and mg are the number density, energy, and mass of the 
target material which are the ground state electrons of the hydrogen atoms, respectively. 
Equation (1) is based on the fact that for hydrogen atoms, the kinetic energy Eg is equal 
to the ionization potential U. 

An appreciable amount of energy can also be transferred to the free electrons that 
are formed as the propellant is ionized. An expression to determine the amoimt of this 
energy exchange is derived in appendix C by use of the kinetic theory of reference 13. 
From equation (C4) 


dx 



In 


^3.45X10® E|^Tgi 


”e“R 


( 2 ) 


where and Tg are the number, density, and temperature of the free electrons, 
respectively. The range of conditions over which the assumption that the free electrons 
thermalize is valid is discussed in appendix D. 

In encovinters of beam particles with neutral target atoms, the energy transfers are 
integrated over a range having the ionization potential of the target material as the lower 
limit and a head-on collision defining the upper limit of integration. The range of impact 
parameters for such encoimters is on the order of the Bohr radius. When considering 
energy transfer to free electrons, however, the minimum energy exchange is cut off at 
an impact parameter corresponding to the Debye length for the case of a charged incident 
particle. The many events in which appreciable energy is exchanged in this relatively 
large realm of influence of charged incident particles make the total energy exchai^e 
from neutral incident particles to free electrons small in comparison. The energy ex- 
changes from neutral incident particles to free electrons will for ease of calculation be 
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(slightly pessimistically) neglected herein. 

Let be the number of species R particles leaving the reactor per unit time. 

The energy transfer rate from reactor species R per xmit distance x to both bound and 
free electrons is thus 
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Conservation Equations of Propellant Species 


Continuity of free electrons . - Consider the control volume in Lagrangian coordi- 
nates given in figure 2. The free electrons are magnetically confined so that, when dif- 
fusion is neglected, the magnetic field lines are envisioned as forming the channel 
boundaries. A coinciding wall of solid material was assumed adequate to prevent loss of 
neutral hydrogen atoms. 

For a steady process, the number of free electrons created by collisions per imit 
time inside the control volume equals the net flow of free electrons out of the control 
volume; that is. 



— (n A Ax) 
6t ® 


coll 


^Ve,xA) 

0X 


Ax 


(4) 


where v „ = v (x) is the mean velocity of the 

6 y X, 6 y 

electrons that had their origin in the propellant. 
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The subscript x denotes that the velocity is in the x direction, that is, along the stream 
tube axis. The free electrons from the beam make a negligible number of ionizing en- 
counters and thus make no contribution to this equation. 

Ionization by beam ions . - The number of ionizing encounters between beam ions and 
hydrogen neutrals per unit time in a distance Ax is 

R 

where Qj^ is the ionization cross section. From equation (C5) it can be written as 


Qr " 


,4 ”^R 
m 


i 

3 12 


2 + ^ for U s KEj^(l - Wjj) 


+ 4)' 


3/2 


(5) 
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Ionization by free elec trons . - The cumulative effect of encounters between free 
electrons and hydrogen atoms can contribute to ionization also. The following expression 
is derived in appendix C for the number of ionizing encoimters per unit volume per imit 
time by free electrons having a Maxwellian velocity distribution: 


VeQ^e = 


n e 
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( 6 ) 


The continuity equation for free electrons can thus be written 

dx 
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Momentum of free electrons . - The divergence of the magnetic field vector 5 for 
cylindrical symmetry is 


lA(rBp.!^=0 

r 9r 9x 

For slow variations of the magnetic field along the axis of symmetry of any local flux 
tube, 



9x dx 


and the first expression reduces to 

2dx 

In the time interval between collisions, the velocity v , of free electrons normal to the 
magnetic field describes arcs of Larmor radii: 


r = 


mv 


e,i 


eB 


( 8 ) 


Substitution of the last two expressions into the Lorentz force equation for the component 
parallel to the B -fie Id yields 
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II 




”^e^e,l dB 
2B dx 


( 9 ) 


The slow B-field change stipulation in this derivation is valid since the electrons make 
many orbits during only a very small B-field variation in the application herein. Equa- 
tion (9) is then valid if the collision frequency is small compared with the cyclotron 
frequency (see appendix D). 

Integration of equation (9) over a Maxwellian distribution in the transverse plane 
yields an average force per electron of 


F|| = 

B dx 
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This force combined with any electrostatic force due to charge separation -q(d(p/dx) 
and the frictional forces due to collisions with hydrogen ions and atoms is balanced by the 
inertial force on the electrons. Thus, 


dv 




e,x _ 
dx 


kT 


e ^ 
dx 


- “e‘'e,n(^e.x ‘ Vx> d”) 


Several underlying assumptions permit the momentum equations to be written in this 
simple form. Collisions between like particles result in a pressure term that was 
negligible for the cases of interest when compared with the large electromagnetic terms. 
Therefore, it is included only in the momentum equation for neutrals. Collisions between 
unlike particles are accounted for by the "friction terms. " These terms are written in 
the phenomenological form of SchlUter (ref. 14), which involves a difference of mean 
velocities of the various species and an average momentum transfer collision frequency 
u. Also the three propellant species are assumed to have nearly equal mean velocities. 
This assumption, which will be borne out in later calculations, eliminates a small inertial 
term accounting for the mass source of a given species due to ionization and charge ex- 
change. 

In fact, the coupling between free electrons and hydrogen ions will be assumed to be 
so close that the velocities (and thus accelerations) of these two species in the x-direction 
are equal (v„ „ = v. = v ). Thus, the first friction term in equation (10) is eliminated. 

0 ) X. X f X X 

This process of accelerating ions by their strong coupling with electrons is sometimes 
termed "adiabatic acceleration" and is discussed in references 15 to 17. This coupling 
is not possible miless the electron paths for recirculation, or short circuiting, are elimi- 
nated. Thus, for example, any wall used for additional confinement must be nonconducting. 

The collision frequency for the remaining friction term was evaluated by the follow- 
ing expression (see pp. 57 and 59 of ref. 13): 


e, n 




where g = |^ 


- Vj^ I , and as discussed on page 10 of reference 18 
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Thus for the study herein 


a 



The momentum transfer cross section a is equal to the total scattering cross 

0 9 n 

section when the "hard sphere” model for neutrals is assumed. A constant value of 
-19 

0. 7X10 square meter was considered sufficiently representative for ^ in the 
energy range of interest (fig. 1. 5 of ref. 19 and p. 199 of ref. 20). Integration of the 
preceding expression then resulted in 


V = 2<10 

e,n 


-12 




( 11 ) 


Momentum of protons . - The hydrogen ions (protons) receive very little direct energy 
from beam particles (see the section entitled Stopping power). They are accelerated by 
the electrons throu^ the electrostatic potential cp and retarded by their collisions with 
neutrals. The adiabatic acceleration assumption eliminates the friction term between 
protons and electrons. The momentum balance in a parallel direction is 


- *^^i X 
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dx 
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where in this case 
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and 



An expression for scattering cross section was obtained by fitting a hyperbola to the data 
in table I of reference 21 (see fig. 3): 
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The corresponding expression for collision frequency is 


p.= 0.867X10“^^ n„ 
i,n n 


(13) 


Eliminating the electrostatic potential term by substituting equation (10) into equa- 
tion (12) and using expressions (11) and (13) for collision frequencies result in 


dv kT ( 

m.V — ^ = - — - — - m.n„(v„ - v„ JIO. 433x10"^'^ + 1.09x10 
^ * dx B dx 1 n X n,x \ 




where m. » m^. 
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Neglecting diffusion across field lines yields the following relation between flow area 
and the magnetic field: 



Thus 


-i ^ ^ 

B dx A dx 

and the preceding momentum eqiiation becomes 

dv kT ^ / 

dx v„Am. dx "\ 
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Momentu m of n eutral propellant atoms . - The momentum balance for neutrals can be 
written 
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Using the previous procedure for evalviating collision frequencies (with o- _ = a 
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n n n 


where the equation of state was applied with cold neutrals to the pressure gradient term. 

Calculations indicated (Preliminary Results section) that the friction force between 
the ions and neutrals is very effective in accelerating the neutrals for cases of interest 
herein, and thus v^^ » v^. The acceleration of neutrals would actually be aided further 
by charge exchange (neglected herein) between these two species of the propellant. 

Propellant energy . - The energy added to the propellant by beam particles minus the 
energy lost by ionization and radiation equals the energy transported away by the propel- 
lant. Thus 
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is the amount of energy to heat the newly formed free electrons to Tg. The loss due to 
bremsstrahlung is (p. 31, ref. 22) 


= 5. 35x10-3’ „^(n^ + „^) 


where the loss is expressed in Newton-meters per cubic meter and where electron tem- 
perature is expressed in Mloelectron volts. 

^^1 il I I I ^ The cost of forming a hydrogen ion U' by 

the free electrons is quite high since their 

■ average energy is even less than the ionization 

potential. In this low energy range the proba- 
bility of exciting a hydrogen atom to a low energy 
level can be much greater than that of an ionizing 
encounter. Calculations based on the cross- 
sectional data of reference 23 and averaged for 
a Maxwell distribution of incident electron veloc- 

^ ity are given in figure 4. The following empiri- 

^ ycu'ateddata from ref. 23^q | ^ | (shown by the dashed line) will be 

0 . 002 . 004 . 006 . 008 . 010 . 012 used in the remaining calculation: 


Maxwell average of cal- / ' 

“culated data from ref. 23^ 

I I I I I 1 
.002 . 004 . 006 . 008 . 010 

Electron temperature, Tg, keV 


Figure 4. - Cost of ion production by low-temperature 
electrons. 
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U* = 0.027211 + 0-0Q34 ] 


For the relatively high energy ion beam, radiative encounters were neglected, and 
U' was set equal to 0. 0136 kiloelectron volt. 


Conservation of Propellant Mass 

The flow rate of the total propellant species m is 

m = [(n^m^ + n.m.)v^ + nnVn,x]^ ("e“i^x + %Vn,x)^ = constant 


(18) 


since n^ = n^ and « m.. 


APPLICATION 


The preceding equations will now be applied to the species and ranges of conditions 
pertinent to a conceptual thermonuclear rocket. 


Reactor 

In the reaction of interest for propulsion applications (refs. 5 and 6) the two isotopes 
helium 3 and deuterium combine to give off a 3. 5 MeV alpha particle and a 14. 7 MeV 
proton. If the energy of these reaction products is expended in heating the incoming 
fuel, then the number of fuel particles inside the reactor should outnumber the reaction 
products by a ratio of about 100 to 1. The beam emitted by the weaker mirror would then 
be essentially composed of only fuel particles. 

The flux and initial energies of the beam are dependent on the magnetic mirror 
system and its particle loss mechanism. The author of reference 24 calculated a large 
electrostatic potential about a reactor, which was caused by the greater tendency of the 
electrons relative to that of the ions to scatter into the "loss cones" of the mirrors. The 
kinetic energy of the electrons is greatly reduced as they emerge, whereas that of the 
escaping ions could be easily doubled by the electrostatic field. This additional kinetic 
energy is in a translational mode, the velocity of which is parallel to the axis of symmetry 
of the reactor. 

Of all the ions inside the reactor, however, the lower energy ones are most apt to 
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TABLE I. - DATA FOR REFERENCE CONFIGURATION 


Ion 

Electron 

Energy of ions 

Flux of ions 

Mass flow 

Minimum 

Initial 

temper- 

temper- 

emitted from 

emitted from 

of beam, 

inside 

cross 

ature 

ature 

weaker mirror, 

weaker 

kg/sec 

diameter 

section 

inside 

inside 

®He-3n 

mirror. 


of mirror 

of ion 

reactor. 

reactor, 


%e-3 ^ 


coils. 

beam 

keV 

keV 

keV 

-1 


^m, min’ 

(jet). 




sec 









m 

^ 0 ’ 







m 

50 

38.5 

50 

0 . 625x10 ^^ 

0 . 522x10 "^ 

1.22 

0 . 0942 

100 

67.0 

100 

.312 

.261 

1.22 

.0805 


enter the loss cone (p. 351, ref. 22). Thus, the translational energy of the escaping ions 
is somewhat near their mean thermal energy inside the reactor. In order to size the 
reactor these two energies will be assumed equal; however, because the precise relation 
is uncertain, the calculations herein will be identified by values of initial beam energy 
instead of reactor temperature. 

The reactor used is essentially that described in reference 5 and was sized to supply 
100 megawatts of power to the escaping ions. Pertinent data are shown in table I for a 
reference configuration. The reaction zone for this reference configuration was arbi- 
trarily given a length -diameter ratio of 5. For the given total power, this ratio deter- 
mines the reactor diameter. A shield was inserted between the reactor and the magnet, 
thus influencing the min imum magnet diameters. A mirror ratio of 4 was used imless 
stated otherwise. 


IVlagnetIc Field Strength 


It was assumed that the reactor particles and propellant do not alter the imposed 
magnetic field. The field near the axis of symmetry of the magnet (fig. 5) was obtained 
from reference 25. An empirical curve (dashed) was used to express B/Bq and A/Aq 
as a fimction of x near the axis of symmetry. The resulting expression is 




-0.231 



+ 2.531 



for 0 i 



(19) 
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Figure 5. - Magnetic field downstream of reactor mirror 
on axis of symmetry. 


Reduction of Equations 

Substituting for the constants in equations (3), (7), (15), (16), (17), and (18), re- 
arranging, and considering two species of ions from the lOO-megawatt reactor result in 
the following set of equations: 
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■14 


+ 0.2512 
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2 \ 3.21 / 
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1 rjrjrj _ 0. 1801 / 74.8 _ 17.30 ^ 


/^He-3 \^He-3 V^He-3 


when 3. 21 < Ejj^ g < 108. 8 


0 when 3 <3.21 


®+ii^ when Ep. a 72.5 
3 ./eT ® 
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+ 
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H" 


f 


1 - + 12. 667 ../l - .+ 2.221 F -2. 


H" 


H" 

E, 


F ^ = F ^ + 0. 375 
D+ 


(l - F A 

V H7 100 


n 




H + n„ 
e n 



(24) 

(25) 

(26) 


The energies of the reactor particles are given in kiloelectron volts instead of joules 
in these final equations. 


RESULTS AND DISCUSSION 
Preliminary Results 


Simultaneous solution of the preceding conservation equations for the propellant and 
the beam was made on a high-speed digital computer. 

The early calculations indicated that v „/v was so close to unity that the com- 

A A 

puting process was made difficult because of the loss of significant figures when subtract- 
ing V /v from 1 in equation (20). Physically this was due to the very effective accel- 
eration of neutrals by the hydrogen ions for the cases of interest. In all later calculations 
v_ „/v„ was set equal to 1. Equations (15a) and (16a) were, therefore, replaced by 

11^ X X 



■>6*^6 dA 
m dx 


Computations were also facilitated by assuming a slightly nonzero initial fractional 


ionization of the propellant F . . Results were insensitive to a range of initial F 


from 10"*^ to 10"® 


H" 




A value of 10 was thus selected. 

To establish an initial condition, it was assumed that the flow of hydrogen propellant 
was choked at the injection station. The results did not change appreciably as the hydro- 
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Energy of deuterium, Ep, keV Velxity of hydrogen, v^, m/sec 


ISD 

CO 
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1 


(e) Energy of deuterium in beam. 


Electron temperature, Tg, eV Energy of helium 3, keV 



B, X, m 

(f) Temperature of propellant electrons, 

Concluded. 


gen temperature leaving the shield was increased from 1000° to 2000° C. A temperature 
of 1000° C was used in the remaining calculations. 

Preliminary calculations indicated that the temperature of the propellant electrons 
was so low that loss of energy by bremsstrahlung was negligible. 


Dependence of Typical Results on Distance x 

A convenient way of varying magnetic field gradient while retaining equation (19) is 
by varying D^. The variation of propellant and beam parameters with distance x is 
shown in figure 6 for four different magnet diameters. Ionization of the propellant is 
most rapid at low values of x where the propellant density is hi^ and where large 
amoimts of energy can be transferred from the beam (fig. 6(a)). Density decreases very 
rapidly with distance x (fig. 6(b)) because of the increasing velocity (fig. 6(c)) and the 
increasing cross-sectional area of the flux tube. Beyond a distance of about one magnet 
diameter downstream of the injection station, the velocity increased, and the target den- 
sity decreased to such an extent that little additional energy transfer took place (figs. 6(d) 
and (e)). The fractional ionization then stays essentially constant; however, accelera- 
tion continues for about two more magnet diameters because of the continuing high gradi- 
ent of the magnetic field (fig. 5) interacting with the remaining energy n kT in the 
propellant electrons. The decrease of electron temperature in this region (as the random 
energy of the electron gas is converted to a rearward directed mean motion) can be seen 
in figure 6(f). 

The large effect of magnet size demonstrated in figure 6 is due to the fact that a 
larger magnet permits a more gradual decrease in field strength and thus a slower de- 
crease of target density with distance x. Therefore, greater amounts of energy are 
transferred from the beam. The smaller field gradient of a larger magnet gives less 
propellant acceleration at low values of x; however, acceleration continues over a 
greater distance. This extended acceleration and the larger amount of energy exchange 
result in a net gain in propellant velocity at large values of x. 

Thrust is defined herein as the component of momentum in the direction of the axis 
of symmetry and can be written as 


•^= <“He-3^He-3 + + “V°°®av^ 


The residual momentum of the reactor particles 


™He-3^He-3 + 
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Infttal flow area. 


1 2 3 

Distance ratio, xlD^ 

Figure 7. - Variation of flow divergence with 
distance. 


Distance ratio, x/D^ 

Figure 8. - Dependence of thrust on the axial 
station at which it is calculated. Magnet 
diameter, 2.44 meters; propellant flow rate, 
0.01 kilogram per second; initial beam 
energy, 50 keV; initial beam area, 

0.0942 square meter. 
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(a) Effect of magnet size. Initial beam area, 
0.0942 square meter; initial beam energy, 
50keV. 
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(b) Effect of initial beam area. 


Figure 9. - Jet performance at hydrogen flow rate of 0.01 kilogram per second. 





is of significance only at the lowest flow rates of hydrogen. The average cosine of the 
angle between the propellant flow and the axis of symmetry was weighted on a basis of 
area normal to the flow. The expression used ‘herein, exact for ideal conical nozzle flow 
(ref, 26), is 


1 + cos 0Vi 

cos e = 

av 2 

where 6-^ is the angle between the outer or bounding fl^IX line and the axis of symmetry. 
Plots of cos^y 6 against distance x are shown in figure 7. Beyond a value for x/D 
of 3, the decrease of cos d quite closely balances the small increase of v„. Thus 
beyond this distance it makes little difference at what x station thrust is calculated 
(fig, 8). Results should therefore be quite insensitive to the location at which the high- 
velocity jet eventually diffuses across the weak field lines and detaches from the magnetic 
field. 


Jet Performance 

Specific impulse is defined in terms of thrust and total mass flow in the jet as 


(niHe-3 + “d 

The product of times the acceleration of gravity is the mean effective velocity of 
the jet. 

The dependence of ^ and > on magnet size and initial beam area Aq is shown in 
figure 9. Figure 9(a) is based on the data of figure 6 at a distance x of 10 meters. The 
large influence of magnet diameter is again demonstrated. The only influence of 
magnet strength in the calculations is through its effect on the size of the cross-sectional 
area of the beam and magnetic flux tube. Increasing the mirror ratio for the same 
reactor diameter decreases Aq according to equation (14). Also for the same reactor 
volume, Aq can be decreased by using a longer and more slender reactor (greater 
length -diameter ratio L/Dj^). For example, Aq can be cut in half by either doubling 
the mirror ratio or by increasing the length -diameter ratio of the reaction zone 
by a factor of 2.8. Various pertinent combinations of Bq/Bj^ and are indicated 

on the abscissa of figure 9(b). Reducing the cross-sectional area of the interaction 
region for the same hydrogen flow, beam flux, and energy results in an increase in energy 
transfer. In general, figure 9(b) shows that decreasing Aq increases both ^ and 
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Hydrogen flow rate, m, kg/sec 

Figure 10. - Effect of propellant flow rate on jet performance. Initial beam energy, 50 keV; magnet 
diameter, 3. 66 meters; initial beam area, 0.0942 square meter; distance, x, 10 meters. 

over the range of variables studied. 

The effect of initial beam energy can be seen by comparing figures 9(a) and (b). The 
four data points in figure 9(b) represent performance at an initial beam energy of 50 keV. 
The beam area for the 100 -keV case would have to be reduced to about one -third of its 
original value to bring the thrust and specific impulse up to that of the 50-keV case. 

These results t37pify the inverse dependence of energy transfer rate on beam energy level 
as determined by equation (3). 

The influence of propellant flow rate on jet performance is shown in figure 10. In- 
crease of m increases ^ but decreases As m is increased, the target density is 
increased to a condition where the beam is able to transfer essentially all its energy to 
the target. Beyond this point the acceleration of the hydrogen continues as the magnetic 
field gradient converts the residual n^kT^ to directed translational energy. The location 
where the beam runs out of energy is apparent by the abrupt change in slope of the curves 
in figure 10. To save computing time, the calculations were completed with the remain- 
ing neglibible energy set equal to zero, when the energies of both species of reactor 
particles fell below 0. 05 keV. 

The power that is useful for propulsion is often defined as . The ratio of 

this jet power to initial power in the beam is shown in figure 10. Also shown is the 
kinetic energy in the hydrogen alone, the energy remaining in the beam that could not be 
transferred to the hydrogen propellant, and the cost of ionizing the propellant, which 
includes radiative encounters. For this case, the maximum overall efficiency of the 
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Figure 11. - Collision parameters. Initial energy of 
beam, 50 keV; initial magnetic field strength, 

20. 8 Webers per square meter; hydrogen flow rate, 
0.010 kilogram per second; magnet diameter, 

2.44 meters. 



Distance, x, m 


Figure 12. - Variation of diffusion parameter with dis- 
tance. Initial energy of beam, 50 keV; magnet diam- 
eter, 2.44 meters; hydrogen flow rate, 0.010 kilogram 
per second. 


energy transfer and acceleration process was calculated to be 25 percent. This maximum 
occurred near an m of 0. 01 kilogram per second. At this point of maximum efficiency, 
the jet thrust was more than seven times that of the beam alone. At low hydrogen flow 
rates only a small amount of power from the beam could be transferred to the target. At 
high flow rates loss is mainly due to power consumed in ionizing the hydrogen. 


Appraisal of Assumptions 

The computed data, such as that presented in figure 6, permit appraisal of several of 
the assumptions regarding collisions discussed in the ANALYSIS section. Details are 
given in appendix D and results are plotted in figures 11 and 12 for a typical set of condi- 
tions. Parameters are presented that are useful for judging the validity of the assump- 
tions of binary collisions, thermalization of free electrons, beam equilibration, and the 
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expression for the magnetic force on free electrons. As discussed in appendix D, results 
are, in general, quite favorable to the methods used in the analysis. Loss of charged 
particles by ambipolar diffusion across the flux tube (fig. 12) appears to be a much lesser 
problem than that occurring inside the reactor. The large values of diffusion coefficient 
for large values of x permit the high-velocity jet to diffuse eventually across the mag- 
netic field lines, which it must do to produce thrust. 

CONCLUSIONS 

The acceleration of a hydrogen propellant along a magnetic tube of flux by use of 
high-energy beams of li^t -weight ions was studied. Subject to the assumptions and 
simplifications in the analysis, and for the range of variables studied, the following re- 
sults were obtained: 

1. Energy transfer was most effective when the beam was concentrated on a dense, 
closely confined target. Thus, increase of propellant mass flow and decrease of initial 
beam cross-sectional area improved results. Ionization of the hydrogen was usually 
quite rapid near the station where the beam first impinged on the target. For the cases 
studied, this was also the location of minimum cross section for the confines of the pro- 
pellant flow. 

2. For the case of a magnetic field generated by a short solenoid, the divergence, or 
spreading, of the flux tube depends to a large extent on the solenoid diameter. Increasing 
the diameter of such a magnet permitted suitable energy exchanges over longer distances 
of beam travel. Beyond a distance of about one magnet diameter downstream of the mini- 
mum cross-sectional station, the velocity increased and density decreased such that 
little further energy transfer could take place. Some acceleration of the hydrogen con- 
tinued for about another two magnet diameters because of the continuing high magnetic 
field gradient interacting with the remaining energy of the propellant electrons. 

3. An optimum hydrogen flow rate was fovind where the momentum (thrust) of the 
propellant could be increased to seven times that initially in the beam and with an overall 
efficiency of energy transfer of 25 percent. At low hydrogen flow rates, only a small 
amount of energy could be transferred from the beam to the target. At high flow rates, 
loss was mainly due to power consumed in ionizing the hydrogen. 

4. Bremsstrahlung was negligible since electron temperature was always below 

11 eV. Sizable radiation losses were calculated, however, because of the high number 
of excitation encounters between these electrons and the hydrogen atom. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 24, 1966. 

120-27-06-04-22. 
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APPENDIX A 


SYMBOLS 


[The mks system of units is used unless otherwise specified J 


A 

cross-sectional area of jet 

Q 

ionization cross section 

B 

magnetic field intensity 

r 

^r /(^ly at a reactor temperature 

% 

specific heat at constant pressure 


/ of 100 keV 


for an ideal electron gas, 
0. 441x10^® J/(kg)(keV) 

T 

temperature, keV, imless other- 
wise stated 

D 

diameter 

t 

time 

m 

magnet diameter (see fig. 5) 

U 

ionization potential for atomic 


coefficient of ambipolar diffusion 


hydrogen, 0. 0136keV 

E 

energy per particle 

U’ 

cost of producing a hydrogen ion 

e 

unit of charge, 1.60206x10"^^ C 

V 

velocity 

F 

equilibrium fraction 

W 

defined by expression (20) 


thrust, N 


ratio of target particle (electron) 
velocity to incident R particle 

f 

flux of particles 


velocity 


acceleration of gravity 

X 

defined by expression (21) 

g 

magnitude of relative velocity be- 
tween two particles 

X 

distance from station of minimum 
jet area (and hydrogen injection 

> 

specific impulse 


station) along a magnetic field 

K 

4memj^/(me + mj^)^ 


line 

k 

Boltzmann constant, 1.6x10"^® 

Y 

defined by expression (22) 


joules/keV 

Z 

defined by expression (23) 

L 

reactor length 

z 

charge 

m 

mass per particle 

a 

mean thermal speed 

m 

mass flow rate (without subscript 
it is for the hydrogen propellant) 

€ 

capacity of a vacuum, 
8. 854X10"^^ F/m 

n 

number density 

e 

angle between a streamline and 

P 

pressure 


the axis of symmetry 
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X ratio of mean to minimum dis- 
tance between particles 

Vg collision frequency between 

jL f 1X1 

species I and m 
momentum transfer cross section 

Jt, m 

between species S. and m 
<p electrostatic potential 

CO cyclotron frequency 

Subscripts: 
av average 

brem bremsstrahlung 
coll collision 

cx charge exchange 

D deuterium 

e electrons 

H hydrogen 

He helium 

He -3 helium 3 


I ionization 

i ion 

max maximum 

min minimum 

n neutral target particle 

R beam or reactor 

t target 

X component in x direction 

0 station where beam first impinges 

on the target, x = 0 

1 perpendicular to magnetic field 

lines 

1 1 parallel to magnetic field lines 
Superscripts: 
o neutral 

total time differentiation 
~ mean mass motion 

— vector quantity 

+ single charge 
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APPENDIX B 


EQUILIBRIUM COMPOSITION OF BEAMS 

The following derivations follow quite closely the theory in section B of reference 9 
except that in the present case a partially ionized target was used. Also the cross- 
sectional area and velocity of the beam herein are not constant. The continuity equation 
for the singly charged helium beam species, for example, can be written as 


^He^ 


dx V He^ He / 


n 
He° 


n 


+ cr. 


H-* 


r^°’ + n 


He^ 


..n 




n . - OTc^°n n 
H+ He+ 




where, for example, is a charge exchange cross section in which the charge of the 

incident particle is changed from 0 to 1. The interacting species to which such a cross 
section corresponds are labeled on the immediately following number densities in this 
equation. The subscript I denotes an ionizing encounter. It was assumed that the inci- 
dent particles undergo frequent enou^ charge exchange and ionization events that 

V = V . Also the propellant velocity was usually so much lower than the beam 
He° He'^ 

velocity that the former was neglected. 

»Q. 

Sufficient data to adequately determine Uj were not available. It is therefore 
assumed herein (as in the ANALYSIS section) that a li^t neutral is as effective in under- 
going an ionizing encounter with a neutral atom as it is after taking on a single charge. 
Thus, in the above equation 


If the beam equilibrates at a certain composition so that the number of neutralizing 
events is balanced by the munber of ionization encoimters, then 

Dividing the right side of the former equation byfn ,+n _)fn ,+n j and defining 

V He+ He°/\ H°/ 


n 


He^ 


= fractional ionization of the helium beam 


He n . + n ^ 
He^ He° 
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and 






s fractional ionization of the hydrogen propellant 




result in 


where 


CT 


0 ,+ 

cx 


F 

He 


o 



+ cr. 


o,+ 


F 

He 


o 


a'*"’ °F F 

He'*' H° 


= 0 


F 

He 


o 



= 1 




= 1 


Solving for F . yields 
He'*' 



This expression approaches the result of reference 9 as F , approaches zero. The 

H'*' 

data of references 9 and 27 were then used to determine the preceding cross sections in 
terms of the energy of the helium beam. The results are shown on figure 13(a). The 
basic parameter is the relative velocity between incident and target particles and not the 
energy, which is more convenient to use. The corresponding energy of helium 3 particles 
is also shown in figure 13(a), wherein it is assumed that the absence of one neutron in 
the helium nucleus does not alter the basic interaction based on relative velocities. The 
calculated data are on a per atom basis while the experimental data are for a target of 
molecular hydrogen. 

Similar results can be obtained for protons in hydrogen by using the data of refer- 
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(a) Helium beams in hydrogen target. (b) Hydrogen beams in hydrogen target. 

Figure 13. - Equiiibrium fractional ionization of beam. 


ence 28. These results were used for a deuterium beam impinging on hydrogen. In these 
cases, however, the charge exchange cross sections should be symmetrical; that is, 

0-°’ which resulted in 



a®’ ■*'F + CT?’ ^ 

CX JJ+ I 


CX I 


(B2) 


(Results are shown in figure 13(b). ) 

The following empirical fits were made to the final results for convenience in later 
computations: 
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( 26 ) 


p . = p . + 0. 375(^1 - P \ ^ 

D"*" H"’' V !!■*■/ 100 

These expressions are shown by the dashed curves in figure 13. The empirical fits 
deviate as much as 15 percent from the calculated results at some of the higher beam 
energies. This, however, is close to the experimental scatter of the basic experimental 
cross sections. 
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APPENDIX C 


ENERGY TRANSFER AND IONIZATION CROSS SECTIONS 

The theory of classical particles is used to determine the energy transfer from beam 
ions to the electrons boxmd in the hydrogen atoms and also to the free electrons liberated 
from the hydrogen atoms. This theory is also used to determine cross sections for ioni- 
zation of hydrogen by the beam ions and by the free electrons. Since all incident parti- 
cles of this appendix have a positive charge, the + superscript used in the ANALYSE 
section has been deleted. 


Energy Transfer to Bound Electrons 

In this case the mass of the target particle m^ is much less than that of the incident 
particle mj^. The cross section ct(AE) for the transfer of the amount of energy AE is 
given by equation (24) of reference 2. After correcting several algebraic errors (see 
appendix of ref. 29) it can be written in mks units as 


cr(|AE|) 


2 4 “r 
”e 


E 

1 + i ® 


3 |ae| 


when AE < KE„ 1 - — 

' »E 


.3/2 


2 


16jre^Ejj(AE)^{l + — 

4/ 


V 


l3 3 IaEI 


i.iMi (Air 


3|AE|/ 6VrKE^Ej^ 


(Cl) 


when 1 ae| a KEj^ll - — 


Absolute value signs have been used to avoid confusion since reference 2 uses -AE as 
the loss of energy from particle R to particle e, which is opposite to the sign conven- 
tion elsewhere in this report. 

Equation (Cl) is then substituted into the following relation for stopping cross section: 


JL?^= /'^^^^“^a(|AE|)|AE| d(|AE|) 

n„ dx •'U ^ ' ' ' 

n 

where for a head-on collision (eq. (15) of ref. 2) 
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(C2) 




After a tedious integration, the stopping power expression becomes 



Substituting the limits of integration and reducing in terms of v^/vp^ with Eg = U for 
the hydrogen atom result in the following equation for singly charged incident particles: 
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Only a small amount of data (figs. 14(a) and (b)) was found for the stopping power of 
species close to those of interest in the desired energy range. Results of equation (C3) 
should be compared first with the partial stopping power data of figure 14(a), which is 
for charged incident particles. Agreement appears satisfactory for present purposes. 

In the beam composition discussion of the ANALYSIS section, the assumption was 
made that the neutral incident particles are as effective as singly charged incident parti- 
cles in transferring energy to the neutral target particles. This assumption appears 
reasonable but slightly optimistic in figure 14(b). No contribution to stopping power due 
to charge exchange was included in the analysis, which is a somewhat pessimistic assump- 
tion. Comparison of the total stopping power (diamond symbols) with the results of equa- 
tion (C3) indicates that the use of equation (C3) for energy loss rate should be reasonable. 


Energy Transfer to Free Electrons 


From pages 114 and 118 of reference 13 


J_^ 


, 2^4 

1 


k(TR - T^) 


AE 


ire^ /2kT^ 2kT\^/^ 

Vj^m^mnl — + 


In 


max 


AE 


min. 


m 


R 


m. 


The maximum energy exchange expression (eq. (C2)) can be reduced further here since 
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Atomic stopping cross section, (l/n^HdE/dx), keV/m^/atom 
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O Partial (ref. 3)forH'^ 

□ Partial (ref. 4) for 

O Total (ref. 4) 

A Partial (ref. 4) for He° 

Eq. (C3) 

T Limits of experimental scatter _ 


8x10’22 



Proton energy, keV 


(a) Partial stopping power of a proton beam 
In molecular hydrogen. 



(b) Helium beams In molecular hydrogen. 

Figure 14. - Comparison of experimental and 
theoretical stopping power cross sections. 


it occurs in the argument of a slowly varying log 
function: 


4m, 


ATT ^ 

I max 


m. 




R 


The minimum energy exchange | AE was deter- 
mined by setting the impact parameter equal to the 
Debye length (p. 82 of ref. 13) so that 

= 3.45X10® 

^Imin 

-2 

Since, at most, is the order of 10 keV, whereas 
Tjj ^ the order of 10 keV, the energy transfer equa- 
tion reduces to the following for singly charged inci- 
dent particles: 


dE 




R 


4 ^ 
m_ 


dx 




.3/2 


ln| 


r3. 45X10® E^Tg^ 
"e4“R 


(C4) 


Ionization by Beam Ions 

Ionization cross section is defined as the cross section for a collision with energy 
loss greater than the ionization potential U. That is. 


Substitution of equations (Cl) and (C2) into this integral expression results in 
38 




39 


Using Eg = U for the hydrogen atoms and writing energy ratios in terms of mass and 
velocity ratios reduce these expressions to 


Qr = 






+ — (— ) for U:SKE„jl 
3 12\yJ R' 


v3/2 




^R, 


5 ^ 2 V 2 

6 3 


48 V. 



^M!l.2-.41 


'R 



for USKE 1 



(C5) 


Ionization by Free Electrons 

The number of ionizing encoimters per unit volume per unit time between the neutral 
hydrogen atoms and the free electrons having a Maxwellian velocity distribution can be 
written as 


"n"eVe = "n 


(C6) 


where 


f(^2) = “e 


m. 


3/2 


i2jrkT„ 


exp - 


^2 

2kT, 


and 


3 2 

d V 2 = ^JTVg dv 2 for spherical symmetry 

Since interactions are between free and boimd electrons, the subscript 2 will be used 
temporarily for the free electrons and subscript 1 will refer to the bound electrons 
where demarcation is needed. There is no confusion in continuing to let n^ be the 
number density of free electrons and n^^ the nvimber density of neutrals, which is herein 
equal to the number density of bound electrons. 
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The preceding integration is to be made over all velocity space; however, the ioniza- 
tion cross section Q is zero for energy exchanges less than the ionization potential U. 
From equations (26) and (27) of reference 2 


^ ®*2 ®“2 V 



when U + Ej ^ E2 



This expression is for the cases wherein the mass of the target particle is equal to that 
of the incident particle. 

Substituting the Vg = and Ej = U and the last three equations into equa- 

tion (C6) results in 



Numerical evaluation of the two integrals gave the value of Q^v„ plotted in figure 15 as 

© © 

a function of T^. An ©mpirical curv© 


Qe'^e 


e 


4 


47t V2ffmgkTg Ue^ 


0. 0446 

eU/kT 



20,95 

eU/kT 


( 6 ) 
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was fitted to the numerical results. Expression (6) is shown by the dashed curve. 

In an ionizing encounter all the energy given up by an incident electron, beyond that 
used for ionization, is given to the liberated electron in a translational mode. The liber- 
ated electron is now part of the free electron gas, and thus the net energy per imit time 
removed from the free electron gas is just that due to ionization, which is Un^^n^Q^Vg. 
Excitation and corresponding line radiation losses can be accounted for by replacing U 
in this final expression by an appropriate cost per ion U'. 
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APPENDIX D 


APPRAISAL OF ASSUMPTIONS REGARDING COLLISIONS 

Numerical solution of the equations of motion over the ranges of interest now permit 
more detailed appraisals of some assumptions regarding collisions. 


Binary Collision Assumption 

Nearly all the equations developed in the preceding section are based upon the assump- 
tion of binary encounters. This approximation can be applied when the particle number 
densities remain below a certain limit. The reciprocal of the cube root of the number 
density can be used as a measure of the mean distance between nei^boring particles. 

The most severe limitation is when Coulomb forces are involved. To obtain a charac- 
teristic length of interaction for Coulomb encovmters, the author of reference 13 (p. 83) 
used the distance 
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It is approximately the distance of closest approach of two mutually repelling particles. 
This distance is within a factor of ten of the radius of a momentum transfer cross section. 
For the binary collision approximation to remain valid the quantity 
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should be large. 

For the range of application herein X is at least the order of magnitude of 10 for 
interactions between free electrons. A plot of a typical case is shown in figure 11 (p. 27). 
For free electrons colliding with bound electrons, the values for X are greater than 
those presented by the curve except in the very small region when 0 ^ x < 0. 01 meter. 
Collisions involving beam particles are not frequent enough to affect the binary collision 
assumption. 
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Thermalization of Free Electrons 


The time for a distribution of electron kinetic energies to approach a Maxwellian 
distribution and the time between large angle deflections are given in reference 30 (p. 78) 
and reference 22 (p. 96). These expressions are within a factor of 1. 4 of the time be- 
tween the momentum transfer collisions of reference 13, which is given as 


e, e 


e,e 
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The distance for thermalization v^/v. „ is plotted in figure 11 for a typical set of condi- 

c X 6, e 

tions and is less than 10" meter over the most important part of the curve. The elec- 
trons should thus have ample time to thermalize. 


Magnetic Force on Electrons 

As discussed in the ANALYSIS section, the expression for the force that the magnetic 
field exerts upon the free electrons is valid if the time that the electron motion is domi- 
nated by the magnetic field is large compared with the time that the motion is collision 
dominated. The usual parameter for such evaluation is the ratio of cyclotron frequency 
to collision frequency. This parameter was most critical for collisions between free 
electrons for the range of variables studied herein. Equation (Dl) was used to determine 

collision frequency t*- _• A typical case is shown in figure 11 where it is seen that the 

0,6 2 

magnitude of the parameter is approximately the order of 10 . 


Equilibrium of Beam 

As an incident particle traverses the hydrogen target it goes through many cycles of 
capturing and losing an electron. The number of these cycles that occur as the energy 
of the particle is degraded from Ej to E 2 is (p. 1142 of ref. 9) 
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Nc = 


where a is the cross section for cycle completion. This expression provides a 
measure of the ability of the beam to equilibrate. The number N is strongly dependent 
on E only; thus, the conclusions of satisfactory equilibration reached in references 4 
and 9, for example, should apply here also. 

The number of cycles per keV of energy degradation Nj,/AE is plotted for a typical 
case in figure 11. This ratio is high in the region of most importance; that is, where 
most of the ionization and acceleration of the propellant takes place. At a value of x 
of 2. 5 meters, however, it falls to a value of 0. 15, and then increases sharply as dE/dn 
approaches zero. The number of cycles per unit distance decreases from 83 at x = 0 to 
less than 1 beyond x = 1. In the lower cycle frequency range the assumption of an 
equilibrium composition is questionable. As a result, the beam may have charge composi- 
tions closer to that estimated for lower values of x. An inaccuracy would occur in the 
energy exchange calculation between the beam and free electrons only. 



Diffusion Across Magnetic Field Lines 

For a partially ionized plasma, the collisional theory coefficient of diffusion perpen- 
dicular to field lines is expressed as a constant times VmT/B^ (ref. 22). The approxi- 
mate expression for the coefficient of ambipolar diffusion across field lines is 


Q , w. 
amb 
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The ratio of ^gLmh propellant species divided by for diffusion inside 

the reactor is plotted in figure 12 for one set of conditions. Since T^ is a constant in 
the propellant and curves of against x have quite the same trends for all cases 

studied (fig. 8(c)), the curve on figure 12 is considered representative of all cases 
herein. In the first two mirror diameters downstream of the mirror coil, the diffusion 
coefficient of propellant across field lines is less than that of reactor species inside the 
magnetic bottle. Beyond a value of x/D of 4, however, the propellant diffusion coeffi- 
cient is more than 100 times that inside the reactor. This is the region, however, where 
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it is desired that the propellant will cross the field lines and continue as a jet in a near 
axial direction. 

The plasma oscillation result of Bohm (p. 462 of ref. 22) indicates that the coefficient 
of diffusion should be proportional to T /B. This would result in an even lower coeffi- 
cent of diffusion in the hydrogen stream compared with that inside the reactor. 


46 



REFERENCES 


1. Gryzinski, Michal: Stopping Power of a Medium for Heavy, Charged Particles. Phys. 

Rev,, vol. 107, no. 6, Sept. 15, 1957, pp. 1471-1475. 

2. Gryzinski, Michal: Classical Theory of Electronic and Ionic Inelastic Collisions. 

Phys. Rev., vol. 115, no. 2, July 15, 1959, pp. 374-383. 

3. Huberman, M, N, : Partial Ionic Stopping Power and the Energy Expended in Electron 

Capture and Loss Collisions of Protons in Hydrogen Gas, n. Phys. Rev. , 
vol. 127, no, 3, Aug. 1, 1962, pp. 799-804. 

4. Cuevas, J. ; Garcia -Munoz, M. ; Torres, P. ; and Allison, S. K. : Partial Atomic 

and Ionic Stopping Powers of Gaseous Hydrogen for Helium and Hydrogen Beams, 
Phys. Rev., vol. 135, no. 2A, July 20, 1964, pp. A335-A345. 

5. Englert, Gerald W. : Study of Thermonuclear Propulsion Using Superconducting 

Magnets. Engineering Aspects of Magnetohydrodynamics, Norman W. Mather and 
George W. Sutton, eds. , Gordon and Breach Science Publ. , 1964, pp. 645-671, 

6. Hilton J. L. ; Luce, J. S. ; and Thompson, A. S. : A Hypothetical Fusion Propulsion 

Rocket Vehicle. J. Spacecraft Rockets, vol, 1, no. 3, May-June 1964, pp. 276-282. 

7. Rose, David J. ; and Clark, Melville, Jr. : Plasmas and Controlled Fusion, MIT 

Press, 1961, p. 421. 

8. Von Ohain, H. ; Keith, B. C, ; and Hunter, R. E. : Survey of Performance Potentiali- 

ties of Non -Airbreathing Propulsion Systems. Vistas in Astronautics, Morton 
Alperin, Marvin Stern and Harold Wooster, eds. , Pergamon Press, 1958, 
pp. 206-213. 

9. Allison, Samuel K. : Experimental Results on Charge -Changing Collisions of Hydro- 

gen and Helium Atoms and Ions at Kinetic Energies above 0. 2 keV. Rev. Mod. 

Phys., vol. 30, no. 4, Oct. 1958, pp. 1137-1168. 

10. Butler, S. T. ; and Buckingham, M. J. : Energy Loss of a Fast Ion in a Plasma. 

Phys. Rev. , vol. 126, no. 1, Apr. 1, 1962, pp. 1-4. 

11. Leighten, Robert B. : Principles of Modern Physics. McGraw-Hill Book Co. , Inc, , 

1959, pp. 494-501. 

12. Evans, Robley D. : The Atomic Nucleus. McGraw Hill Book Co. , Inc. , 1955, 

pp. 575-580. 

13. Burgers, J. M. : Selected Topics from the Theory of Gas Flow at High Temperatures 

(V). The Application of Transfer Equations to the Calculation of Diffusion, Heat 
Conduction, Viscosity, and Electric Conductivity. Parts I and H. Tech. Note 
BN-124a and b (AFOSR TN 58-427 and a). University of Maryland, May 1958. 


47 


II I I 


III 


14. Cowling, Thomas G. : Magnetohydrodynamics. Interscience Publishers, 1957, 

pp. 99-112. 

15. Hendel, H. W. ; and Reboul, T. T. : Adiabatic Acceleration of Ions by Electrons. 

Phys. Fluids, vol. 5, no. 3, Mar. 1962, pp. 360-362. 

16. Rose, D. J. : Acceleration of a Neutralized Ion Beam. Quarterly Progr. Rep. 

No. 53, Massachusetts Inst, of Tech. , Electronics Research Lab. , Apr. 15, 1959, 
pp. 3-4. 

17. Domitz, Stanley: Experimental Evaluation of a Direct-Current Low-Pressure Plasma 

Source. NASA TN D-1659, 1963. 

18. Lyman, Frederic A. ; Goldstein, Arthur W. ; and Heighway, JohnE.: Effect of 

Seeding and Ion Slip on Electron Heating in a Magnetohydrodynamic Generator. 

NASA TND-2118, 1964. 

19. Brown, Sanborn C. : Basic Data of Plasma Physics. MIT Press, 1959. 

20. Barnett, C. F. ; Ray, J. A. ; and Thompson, J. C. : Atomic and Molecular Collision 

Cross Sections of Interest in Controlled Thermonuclear Research. Rep. No. 
ORNL-3113 (Rev.), Oak Ridge National Labs. , Aug. 1964. 

21. Simons, J. H. ; Fontana, C. M. ; Muschlitz, E. E., Jr.; and Jackson, S. R. : The 

Low Velocity Scattering of H'*’ and Hg in Hydrogen. J. Chem. Phys., vol. 11, 
no. 7, July 1943, pp. 307-312. 

22. Glasstone, Samuel; and Lovberg, Ralph H. : Controlled Thermonuclear Reactions. 

D. Van Nostrand Co. , 1960. 

23. Monnin, Carl F. ; and Prok, George M. : Comparison of GryzinsM’s and Born’s Ap- 

proximations for Inelastic Scattering in Atomic Hydrogen. NASA TN D-2903, 1965. 

24. Ben Daniel, D. J. : Plasma Potential in a Magnetic Mirror System. J. Nucl. Energy, 

Pt. C, Plasma Phys. , vol. 3, no. 4, 1961, pp. 235-241. 

25. Callaghan, Edmund E. ; and Maslen, Stephen H. : The Magnetic Field of a Finite 

Solenoid. NASA TN D-465, 1960. 

26. Engle rt, Gerald W. ; and Kochendorfer, Fred D. : Estimated Performance of Radial- 

Flow Exit Nozzles for Air in Chemical Equilibrium. NASA Memol-5-59E, 1959. 

27. Gilbody, H. B. ; and Hasted, J. B. ; Ionization by Positive Ions. Roy. Soc. Proc. , 

ser. A, vol. 240, no. 1222, June 11, 1957, pp. 382-395. 

28. Fite, Wade L. ; Stebbings, R. F. ; Hummer, David G. ; and Brackmann, R. T. : 

Ionization and Charge Transfer in Proton -Hydrogen Atom Collisions. Phys. Rev., 
vol. 119, no. 2, July 15, 1960, pp. 663-668. 




48 



29. Bauer, Ernest; and Bartky, Charlotte E.; Calculation of Inelastic Electron -Atom 

and Electron-Molecule Collision Cross Sections by Classical Methods. Publ. 
No. U-2943, Philco Corp. Aeronutronic Div., Jan. 1965. (Available from DDC 
as AD627033. ) 

30. Spitzer, Lyman: Physics of Fully Ionized Gases. Interscience Publishers, Inc. , 

1956. 


NASA- Langley, 1966 E-2847 


49 


”The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenome?ia in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.'' 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washinston^ D.C. 20546 



